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Abstract 

In this paper, we mainly study a hydro-dynamic system modeling the 
flow of nematic liquid crystals. In three dimensions, we first establish 
local well-posedness of the initial-boundary value problem of the system. 
Then, we prove the existence of global strong solution to the system with 
small initial-boundary condition. 
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1 Introduction 

In this paper we consider the following hydrodynamic system modeling the flow 
of liquid crystal materials in three dimensions (see [TJ [U [3J [5] ) 

|+u-V«- + VP = -AV • (Vd © Vd), in R+ x Q, 

(1.1) { |f +u-Vd = j(Ad+ | Vd\ 2 d), m R+ x SI, 
V • u = 0, in K+ x fi, 

with initial-boundary conditions: 

(1.2) (u(0,a;),d(0,a!)) = (uo(a;),do(a;)) i£(l, 

(1.3) (u(t, x), d(t, x)) = (0, d (x)) (t, x) G R + x fl 

Suppose that f2 C R 3 is a bounded smooth domain, u(t, x) : R + x SI — > M. 3 stands 
for the velocity field of the flow, d(t,x) : R + x Q — > S 2 , the unit sphere in R 3 , 
is a unit-vector field that represents the macroscopic molecular orientation of 
the liquid crystal material and P(t, x) : R + x SI — > R is the pressure function. 
The constants v, A, 7 are positive constants that stand for the viscosity, the 
competition between kinetic energy and potential energy, and microscopic elastic 
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relaxation time for the molecular orientation field. Vd Vc? denotes the 3x3 
matrix whose the entry is given by Vjd • Vjd for 1 < i, j ' < 3. It is easy 
to see that VdQ Vg? = (Vrf) T Vo?, where (Vd) T denotes the transpose of matrix 
Vd. 

System (1.1) is a simplified version of the Ericksen-Leslie model. General 
Ericksen-Leslie model reduces the Ossen-Frank model in the static case, for the 
hydrodynamics of nematic liquid crystals developed during the period from 1958 
to 1968 [2 HI II]- Since the general Ericksen-Leslie system is very complicated, 
we only study a simplified model of the Ericksen-Leslie system which can derive 
without destroying the basic structure. It is a macroscopic continuum descrip- 
tion of the time evolution of the materials under the influence of both the flow 
field u(t, x), and the macroscopic description of the microscopic orientation con- 
figurations d(t,x) of rod-like liquid crystals. The system (1.1)-(1.3) is a system 
of the Navier-Stokes equation coupled with the harmonic map flows. 

In a series of papers, Lin [5] and Lin-Liu [31 E] initiated the mathematical 
analysis of the system (1.1)-(1.3). Since the Erichsen-Leslie system (1.1)-(1.3) 
with | d \— 1 is complicated, Lin and Liu [3J|^ proposed to consider an approx- 
imation model of Ericksen-Leslie system by Ginzburg-Landau functional. More 
precisely, they replaced the Dirichlet functionals 

- / Vd I 2 dx 
2 in 

for d : ft — » S"^ 1 by the Ginzburg-Landau functionals 

for d : 57 — >• R n (e > 0). In 3 , Lin and Liu proved the global existence of 
solutions in dimensions two or three. In [6] , Lin and Liu proved partial regularity 
of weak solutions in dimension three. Furthermore, Lin and Liu in [7] proved 
existence of solutions for the general Ericksen-Leslie system and also analyzed 
the limits of weak solutions as e — > 0. In [§], Hu and Wang give the existence of 
global strong solution and prove that all the weak solutions constructed in [3] 
must be equal to the unique strong solution. 

Recently, Lin, Lin and Wang [5] studied the system (1.1)-(1.3) in two di- 
mensions. They established the global existence and partial regularity of the 
global weak solution and performed the blow-up analysis at each singular time. 
Hong |T3] proved the global existence of the system (1.1)-(1.3) in two dimensions 
independently. 

The aim of this paper is to establish the short-time solution for general 
initial-boundary condition and the global existence for small initial-boundary 
conditions for the system (1.1)-(1.3). 

Notations In this paper, we denote W m ' q (Q) the set of function in L q (VL) 
whose derivatives up to order m belong to L q (fl). For T > and a function 
space X , we denote by L p (0, T; X) the set of Bochner measurable X — value time 
dependent functions f such that t — y\\ f \\x belong to L p (0,T). The function 
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u G V 2 ' (Qt) is defined as 



ueC({o,T],L 2 (n))nw^(Q T ), 

sup || u(t, •) || L 2 (0) + || Vu \\l 2 (Q t )< oo, 
[0,T] 

1-I,P 

where Qt — (0,T) x 0. TTie space D A F represents some fraction domain 

of Stokes operator in L q (see Sect. 2. 3 in flOf). Roughly, the vector-fields of 
i— i p 

D A p ' are vectors which have 2 — - derivatives in L q (£i), are divergence-free, 
and vanish on dfl. B s q (pi) represents the Besov space f!2f which can be regarded 
as the interpolation space between L q pi) and W s+e ' q (tt). From Proposition 2.5 
in 11 Of . we can get 

(i.4) d 1 ^* sX~* 3 (n) n L«(n). 

Since the space variables are in fl, if there is no ambiguity, we write L q (fl), 
W m ' q {Vt), B s qp {tt) as L q , W m > q , B s q p respectively. 

Definition 1.1. For T > and 1 < p,q < oo, we denote by E T ' q the set of 
triplets (u, d, P) such that 

u G C(0,T;D A ~^ P ) n L p (0,T;W 2 ' q {n) n\Yv' q {n)), 
d t u e L P (0,T; L q (n)),V ■ u = 0, 

d G C(0, T; bJ£~' } ) n LP(0, T; W 2 ' q (Q)), d t d G i p (0, T; L*(fi)), 

P G LP(0, T; W 1,9 (fi)), j n Pdx = 0. 

The corresponding norm is denoted by || • H^p-i. 

II {u, d, P) \\e%*= sup || u || !_!,„+ sup || d || 2(1 _ i, + || u \\lp(o.t-,w^i) 

[0,T] D A , P [0,T] B,,„ p 

+ II ^t" \\lp(0,T-Li) + II d IIlp(0,T;IV 2 '9) + II dtd ||lp(0,T;.L«) ■ 



Our main results can be stated as follows. 

Theorem 1.1. Let Q be a smooth bounded domain in R 3 and (1 — ~) • q > 3. 

If u G D \~J' P and d G B^p"^ f~l C 2 < Q (W), tfien 



(1) There exists a Tq > 0, suc/i i/iai i/ie system (1.1) with the initial-boundary 
condition (1.2)-(1.3) has a unique local strong solution (u,d,P) G E T ' ( 
(0,Tq) x 57. Moreover the solution continuously depends on initial data 



in 





(2) For any given unit vector e G S 2 , there exists a S > 0, such that, if the 
initial data satisfies <io|gn = e and 

|| w || + || do - e || 2(1 _i) < 6, 

D A ," B„ iP " nC^°(9f!) 
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then the system (l.l)-(l.S) has a unique global strong solution (u,d,P) G E^ q 
in (0, T) x n for all T > 0. 

We obtain our results in the spirit of [TUJ [5] . The main difficulty is the low 
intcgrability of nolinear item V • (Vd Vd). To overcome this problem, we add 
a certain condition to the initial data such that we can get the estimate of the 
|| Vd |l°°(o, t-l°°)- The paper is written as follows. In Section 2, we give some 
useful lemmas. In Section 3, we prove the local well-posedness. In Section 4, 
we prove the global existence. 

Remark 1.1 In this paper we only prove the results in three dimensions, but 
we point out that our method can deal with the system (l.l)-(l.S) in higher 
dimensions. 



2 Preliminaries 

In this section, we give some useful lemmas which will be used in the sequel. 

Lemma 2.1 QI] Given 1< p,q < oo, u e K^p ^ and f e L p (0,T;L q ). Then 
the Cauchy problem 

Ou 

(2.1) —-Au = f, u\ t =o=u , 

has a unique solution u satisfying 

(2-2) || U ||m/ 1 >p(0,T;L<!) + II u \\lp{0,T;W 2 -i) 

<Cl(ll / ||ZJ>(0,T;L«) + II "0 || 2(1-1)), 

where C\ is independent of uq, /, and T, Moreover, there exists a positive 
constant Ci independent of / and T such that 

(2.3) SUp | U || 2(i-i)< C 2 (\\ f \\lp(0,T;L«) + || U || 2(1-1))- 
te(0,T) B,,p P B,, P " 

Lemma 2.2 [TU] Let VI be a C 2+e bounded domain in R w and 1 < q,p < oo. 
Assume that uo S D Aq v and f G L p (R + ,L q ). Then the system 

f f -Ati + VP = /, f Q Pdx = 0, 

(2.4) <^ V-« = 0, u|an = 0, 

has a unique solution u, P satisfying the following inequality for all T > 0: 

(2.5) |KT)|| , +( ( T || (VP, V 2 u,9tu) ||^, dt)p 

D Aq P Jo 

<C 3 (\\u \\ 1 _i, p +( f || f(t)\\h)h 

D Aq P Jo 

with C 3 = C(q,p,N,n). 

Using Lemma 2.1, we can prove there is a similar conclusion for initial-boundary 
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problem. 

Theorem 2.1 Let ft be a bounded smooth domain in R 3 and (1 — |) • q > 3. 

If u 6 Bg^^ fl C 2 ' a (9J]) and f G U>(0,T;L"), then the initial-boundary 
problem 



(2-6) — - Au = /, u\o pQt = u 

has a unique solution u satisfying 

( 2 -7) || U || 2(1 _i, + || U \\li>(0,T;W 2 .i) + II 9 t U \\ L v( ,T;Li) 

L°°(0,T;B q , p p ) 

<C(|| / \\ L p(Q.T-M) + II "0 || 2(1-1) +Tp II M \\c*-<*(d(i)), 
where d p Q T = (0, T) x <9£1 U {0} x 0. 

Proof. Since uo S C 2 ' a (<9f2), by the standard elliptic theory we get that there 
exists a unique solution u 1 G C 2 ' a (f2) satisfies 

(2.8) Au^O, u 1 \ dn = u . 

2(1— i) 

It is clear that u a — u\ e _Bq iP p and (u — iti)|ao = 0. Using Lemma 2.1, the 
equation 

5u a * I 

^ - A« = /, u| t= o = uo - Ul, 

2(l-±) 

has a solution u 2 £ 5<j, P p and U2(t, x)\gQ — 0. A direct computation shows 
that Hi + U2 is a solution of the initial-boundary problem (2.6). Using lemma 
2.1 and Schauder's estimate, we deduce (2.7). □ 

Theorem 2.2 Let(l-j)-q > 3 andu e L°°(0, T; L>^/' P ). If d G i p (0, T; !U 2 <«)n 

VF 1,p (0,T;L 9 )nL oo (0,T;i? gi p p ) is a solution of the following nonlinear parabolic 
problem 

(2.9) — - Ad- | Vd | 2 d + u- Vd = 0, m (0,T)xQ, 
(2-10) d|a PQT = do, 

w/iere d : fl -> S 2 . T/ien, d |= 1 in [0,T) x O. 

Proof. Multiplying (2.9) by d, we get 

(2.11) 

9(1 ^ ~ 1} - A(| d | 2 -1) - 2 | Vd | 2 (| d | 2 -1) - u ■ V(| d | 2 -1) = 0. 
By the assumption (1 — |) • q > 3 and (1.4), we have 

(2.12) || Vd|Uoo (0 ,T;L~)<|| Vd|| ..a <||d|| 2(1 _l, 

L^(Q,T:B q ^) L(0,T;B q , p p ) 

(2.13) || U ||loo (0iT;I ,oo)<|| U II 2(1 i <|| U || . 

L°°(0,T;B q , p >) L~>(0,T;D Aq *> ) 
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Noticing that (1 — |) • q > 3, we get 

II v(| d | 2 -l) \\m QT) <c || d \\ l °o {Qt) \\ vd || l2(Qt) 

<C II d II 2(i-i) II d \\lp(o,T:W 2 -i) ■ 

L°°(0,T;B q , p p ) 

Since <9 t d G 1^(0, T;L«) and d G L°°(0, T; L°°), we can get that 

(|d| 2 -l)eC([0,T],L 2 ). 

This yields (| d | 2 -1) £ F 2 1,0 (Qt)- Thus, we have that the function | d | 2 -1 
satisfies the following equation 

(2.14) ^ - A/ - 2 | Vd | 2 / - u • V/ = 0, 

(2.15) /|3P Qt = 0, 

in V^'°(Qt)- The uniqueness implies | d | 2 — 1 = in Qt- This proves the 
lemma. □ 



3 Local well-posedness 

In this section, we prove the local well-posedness of the system (1.1) with the 
initial boundary value (1.2)-(1.3). 

Noticing that VP = V(P — J n Pdx), we can assume that 



L 



Pdx = 0. 

n 

Since the exact values of v, A, 7 do not play a role, we henceforth assume 

v = A = 7 = 1. 

Thus, we can rewrite the system (1.1) as 

^ + u ■ V11 - Ait + VP = -V • (Vd© Vd), in R+ x O, 
;; ' -u ■ Vd = (Ad+ I Vd | 2 d), m M+ x 0, 



(3.1) 



9t 



V • u = 0, in 1+xfl, 

J n Pd 2; = 0, 



with initial-boundary conditions 

/ (u(0,a:),d(0,a;)) = (u (a;), d (x)) x G 0, 
1 J \ (u(i, x), d(t, a)) - (0, d (x)), (t, x) G R+ x an. 

Now, we prove the local existence. 

Firstly, we linearize the system (3.1)-(3.2) and construct approximate solu- 
tions. Set (u°(t, x),d°(t, x)) = (u ,d ). Then given (u n ,d n ,P n ) as the solution 
of 
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(3.3) 



^ - Au n + VP" = -n"- 1 • Vu"- 1 - V • (Vrf™- 1 © Vd"" 1 ), m R+ x fi, 



ad 



— - Ad" = -u"" 1 • Vd""^ | Vd"" 1 | 2 d""\ m R+ x 0, 

V • w" = 0, m M+ x 0, 



with initial-boundary data 
(3.4) 



(u"(0,a:),d"(0,a:)) = (u (a:), d (x)) x e SI, 
(u n (t, x),d n (t, x)) = (0, do(x)) (t, x) e M+ x <9Q. 

From Theorem 2.1 and Lemma 2.2, we can obtain that the sequence {(u, n , d", P")} n eN 
belong to E^' q for any T > 0. Moreover we have the following estimate between 
(u",d",P") and (u"- 1 /" 1 ,?"- 1 ). 



(3.5) || u" +1 || +( / || (VP"+\ V 2 U " +1 ,a tU " +1 ) dt)* 

<C(\\ u || „ +( f || u" • Vu" + V • (Vd" Vd") Hi, 

(3.6) 

II II 2(1 -1, + II IliPCCTiW..) + II 9 t d n+1 \\ LP{0 ,T;L*) 

L°°(Q,T;B q , p p ) 

<C(|| do II 2(1 -i, +( / T II -u n • Vd"+ | Vd" | 2 d" dt)5 + Tp II d ||c^(an)). 

B q,p " JO 

Note that w" +1 |ao = 0. Then the divergence theorem implies 

/ Vw" +1 d:r = 0. 
Jn 

Thus, using Poincare's inequality, we get 

II u n+1 \\ W 2, q < C || V 2 u" +1 || L , . 
Therefore (3.5) can be written as 
(3-7) 

II u n+1 || +( f || (VP n+ \u n +\ V 2 u n +\d t u n+1 ) \f Lq dt)* 

L°°(0,T:D Aq p ) Jo 

<C(\\ u || ,_ Kp +( / T || u" • V«" + V • (Vd" Vd") \\* 9 dt)*). 

D Aq P Jo 

Secondly, we present a uniform estimate for the sequence {(u n , d™, P")} ng N. 
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Define 

(3.8) F n (t) = || u n || 2(i-i) + || u n \\Lp(o.t-.w 2 .i) 

+ II d t U n \\lp(0,P,Li) + II VP" ||z,P(0,t;L9), 

(3.9) £7„(i) = || d" || 2(i-i) + II d n \\Lp(o,t-w 2 -i) + II d t d n \\LP(o,t-,L"), 

L°°(0,t;B v , p p ) 

(3.10) H n (t)=F n (t) + E n {t), 

(3.11) F = ||uo|| ( i_i), p , £^ =|l * II 2(i-i,, H = F + E . 

D A q B qp 

1 — - v 2(1 — -) 

Lemma 3.1 Let(l-^)-q > 3 andu (x) E D Aq " , d (z) e P 9 ,p p T)C 2 ' a (dfl). 
Then, there exists a positive To such that the sequence {(u n , d n , P")} n eN is 
uniformly bounded in E™ ■ 

Proof. Using the fact 

Vd I 2 

(3.12) V ■ (Vd Vd) =V( 1 1 ) + Ad ■ Vd 

=V 2 d -Vd + Ad-Vd 

and (2.12)-(2.13), we have that 

(3.13) F n+1 (t) <C(F + || u n - Vu n \\ 

+ || V | Vd n | 2 || LP( o,t;L,) + || Ad n ■ Vd n \\ L P {0 ,fM)) 

<C(F + i? || U n 1^00(0^.^00) || Vu n || L « (0 ,t;i=c ) 
+ II Vd ||l,oo(0,t;Z,=o)|| Ad n HiP^tjL*)) 



<C(F + t» 



.n 1 1 2 



u 



L°°{0,t;D 



+ II d U || 2(1-1) || d" \\LP(0,t;W 2 .i)) 

L°°(0,t;B gtP v ) 

<C(F + t^F 2 (t) + E 2 (t)), 
due to (3.7). Similarly, we obtain from (3.6) that 
(3.14) E n+1 (t) <C(\\ u n -Vd n || III Vd" | 2 d" \\ L P(0,t;L«)) 

<C(E + t? || U n || L oo(o, t;i oo) || Vd" |Uoo (0it;i oo) 

+ t» || Vd" |||oo(0, t;i oo)|| d" \\ L ~ {Q , t ;L~) +t? || d Hc^OQ)) 

<C(P +t»F„(t)f; n (t) + t^; 3 l (t) +tp || d || C a.«(en)). 
Combining (3.13) with (3.14) yields that 
(3.15) 

Hn+i(t) < C{ff + *» + F n (t)E n {t) + El(t)+ || d Hc^OQ)] + E 2 n (t)}. 

Plugging (3.14) into (3.15), we get that 
(3.16) 

H n+1 (t) <C{H +t 1 p(F 2 (t) + F n (t)E n (t)+El(t)+ \\ d \\c.<-(on)) 

+C 2 (E Q + t 1 pF n _ 1 {t)E n _ l (t)+t 1 pEl_ 1 {t)+t L p || d || C 2, Q(ao) ) 2 }. 
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By a direct computation, we obtain 
(3.17) 

H n+1 (t) < C{H + C 2 H 2 + (C 2 + || d Q || c ».<.(fln))t'(fl2-i + H Li + K-i) 

"(9!i)) 

+tp || d ||c 2 . "(do) • 



+t*{Hl + Hl) + (2C 2 H + || d Hc'-cen))*'^-! + #n-i)}> 



where the constant C only depends on SI. Since — >■ H as i — > 0, we can 
assume that there exists a T > such that H n (t) < 2CKHo and H n _\(t) < 
2CKH on [0, T] for some fixed n, where X = 1 + C 2 i? - Then, 

(3.18) 

H n+1 (t) <C{(H + C 2 H 2 ) + C 3 ti(2 e K 6 H% + 2 5 K 5 H% + 2 4 K i H 4 ) 
+ tUc + 2C 2 Ho+ || do \\ C 2, a{dn) )(2 3 K 3 H 3 + 2 2 K 2 H 2 ) 
+ II d \\c^on) ti(2 6 K 6 H% + 2 5 K 5 H* + 2 4 jY 4 /f 4 )}, 
=CKHo + C{C 3 tp (2 6 K 6 H$ + 2 5 K 5 H^ + 2 4 K 4 H%) 

+ t L p(C + 2C 2 Ho+ || do \\ c ^ { on))(2 3 K 3 H 3 + 2 2 K 2 H 2 ) + ti 
+ II d \\c^(on) ti(2 6 K 6 H$ + 2 5 K 5 H% + 2 i K i H 4 )} 

+ t* || do ||c 2 '°(9Q) 

If we choose T <T such that 

(3.19) (C 3 + || d \\ C 2. a(dn) )Ti(2 6 K 6 H^ + 2 5 K 5 H 4 + 2 i K i H 3 ) 

+ T»{C + 2C 2 H Q + || do \\ C 2, a(dn) )(2 3 K 3 H 2 + 2 2 K 2 H ) 
<K - 1, 

and 

i_ 

T o II \\c 2 ."(an)< CHo, 

then 

(3.20) H n+1 (T ) <2CKH Q . 
Since H n+ \{t) is increasing, we get 

(3.21) H n+1 (t)<2CKH , on [0,T ]. 

Arguing by induction, we deduce that {(u n , d n , P n )} n eN is uniformly bounded 
in E% p . □ 

Jo 

Next, we establish the convergence of this approximate solutions sequence 
{(u n ,d n ,P n )} neN . 

Lemma 3.2 Let(l-%)-q >3 7 u £ D 1 A J' P and do G B^^ nC 2 ' a {dfl). Given 
solution sequence {(u n ,d n ,P n )} constructed in Lemma 3.1. Then, there exists 
a positive constant T\ < To such that {(u n , d n , _P n )} n ^^ converges in E%f. 
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Proof. Let 

(3.22) D{u n ) = u n+1 - u n ; D(d n ) = d n+1 - d n ; D{P n ) = P n+1 - P n . 



Define 
(3.23) 

(3.24) 



DF n (t) = || D(u n ) || + || D(u n ) 11^(0,*,^..) 

L°°(0,t:D Aq p ) 
+ || d t D(u n ) \\ LP{0 ,fM) + || VD(P n ) || L , (0l t ii4 ), 

DE n (t) = || £>(<T) || 2(1 _i, + || D(d n ) 11^(0,*;^..) 

L°°(0,t;B„ P v ) 



+ || |U P(0 ,^), 
(3.25) I>ff„(i) =DF n (t) + DE n (t). 

A direct computation shows that the triplet {D(u n ) 1 D(d n ), D(P n )) satisfies 
' ^p- - AD(u n ) + VD(P n ) 

+V( |v<f r 1|!> ) - Ad™ • Vd" + Ad"- 1 • Vd"" 1 , in M.+ x ft, 



(3.26) 



- AD(d n ) 
= -u n • Vd™ In"- 1 • Vd"" 1 
+ | Vd" | 2 d™- | Vd™- 1 | 2 d"" 1 , 



V • D(u n ) = 0, 

I f n D{P»)dx = 0, 
with initial-boundary value 



in R + x ft, 
in R+ x 0, 



(3.27) 
(3.28) 



D(u n )\ t=0 = D(u n )\ dn = 0, 
D(d n )\ t=0 = D(d n )\ an = 0. 
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Using (2.12)-(2.13), we get 

(3.29) || -u n ■ Vu n + u 11 - 1 • Vk"- 1 \\ LP (o.t:M) 

< || u n ■ V( J D( U "- 1 )) \\ L , m ^) + II Vw"- 1 • Div! 1 - 1 ) \\ LPmL «) 

<t' II «" II i-i.p II ^K" 1 ) II 

L~(0,t;2J A / ) £~(o,t;i> A / ) 

+ i"ll^- 1 H , i-i,, II I>(« n_1 ) II / 1 - i . r , 

L°°(0,t;D Aq p ) L°°(0,t;D A(j p ) 

(3.30) || Ad™ ■ Vd" - Ad"" 1 • W™- 1 ||£P( 0l t;£«) 

<IKII , II ^(d™- 1 ) IU*>(<w) 

+ II d"" 1 ||LP ( o, ti i^.«)|| D{d n ^ 1 ) || 2(1 _ 1} , 

L~(0,t;B„ >p p ) 

(3.31) || V 2 d" • Vd" - V 2 ^- 1 • Vd"- 1 11^(0,*;^) 

< II II 2(1-1) || D(d n ^ 1 ) \\Lp(o,t-,w 2 .i) 

L~(0,t;B 5 , p p ) 

+ II cr- 1 ||LP(o, ti ^..) II £>(rf" _1 ) || 2(1 _i) , 

L~(0,t;B„ jP p ) 

(3.32) || tt " • Vd" - w"- 1 • Vd"- 1 ^(ctji.) 
<^IKII , i-i-pll^d"" 1 )!! , 

+ t'||d n - 1 || 2(i-i) || D(u n_1 ) || 

L°°(0,t;B 9 , p p ) L~(0,t;D A(! P ) 

(3.33) || | Vd" | 2 d»- | Vd"- 1 | 2 d"- 1 |Up ( o, t; L«) 

<2*' II d" II 2 2( i-i, II D^- 1 ) || 2(1 _i, 

L«(0,t;S,, p p ) L~(0,t;B 9iP p ) 

+ *' IK" 1 !! , 2(1 -1, II ^(d"" 1 ) || 2( i-l, . 

i0O (°.*i B «.P " ) L~(0,t;B, B p - " ) 

Using (3.21) and (3.29)-(3.33), we obtain that 

(3.34) DE n {t) < Ct* (DE n -i(t) + DF n -i(t)), 

(3.35) DF n (t) < C(*p£)F n _i(t) + D£?„_i(t)). 

Thus, we get 

(3.36) £/d„+i < 4C 2 (ti +tp)ff n _ 1 (t), 

where the constant C only depends on uq, do and Q,. Choose T\ < T such that 
AC 2 {T{ +T{)<\. Then 

(3.37) DH n+1 {t) < ±DH n -!(t), on [0,Ti]. 

It is clear that {(«", d™, P")}„ e N converges in . □ 

Let (it, d,p) be the limit of {(u n , d", P")}, we prove that (u, d, P) is a solution 
to the system (3.1)-(3.2). It suffices to show that the items in the system (3.3) 
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converge to corresponding items in L p (0, T\ \ L q ). We only prove the convergence 
of V • (Vd" Vd"), the others can be proved similarly. 

V • (Vd" Vd") - V • (Vd © Vd) ||lp(o,t i; l«) 
I Vd" I 2 



(3.38) 



< 



o, 



2 2 

+ || Ad" ■ Vd" - Ad ■ Vd IU^o.t^w) 
<2(|| d" - d ||lp(o,t i; iv 2 .9) II d" || 2(i-i) 

L°°(0,Ti;B 9iJ> p ) 

+ II rf I|lp(o,t i; w 2 .9)II d™ - d || 2(1 _i) ). 

L°°(0,T i; B,, p p ) 

Using the convergence of («",d") and (3.21) we obtain 

(3.39) || V • (Vd" Vd") - V • (Vd Vd) || Lp( o,t i; l ? )- 

(3.40) || u n • Vw" - uVu \\lp(o,tv,l«)^ 0, 

(3.41) || u n • Vd" - u ■ Vd |Up( 0> t i; l<o^ 0, 

(3.42) HI Vd" | 2 d"- | Vd | 2 d ||£p(o,t i; i;«)-> 0. 

It is easy to see Vu = and J Q Pdx = 0. As a direct consequence of Theorem 
2.2, we get | d |= 1. Thus, we have proved that the function triplet (u, d, P) is 
a solution to the system (3.1)-(3.2). 

Now, we prove the uniqueness. Let (u 1 ,d 1 ,P 1 ) and (u 2 ,d 2 ,P 2 ) be two 
solutions of the system (3.1) with the initial-boundary conditions (3.2). Denote 



Su = 



Sd = d 1 -d 2 , SP = P 1 -P 2 . 



Then, the triplet (Su, Sd, 5P) satisfies the following system 

ASu + VSP 
-u 2 - V^ + u 1 ■ Vu 1 - V(^f^) 



(3.43) 



dSu 
dt 



+V(^J_)-Ad 2 -Vd 2 + Ad 1 -Ad\ in (o } Ti)xfi, 



+ | Vd 2 | 2 d 2 - 
V • Su = 0, 



Vd 1 | 2 d 1 



in (0,Ti) x 
m (0,Ti) x 0, 



with initial-boundary value 

(3.44) Su\ t=0 = Su\ 9n = 0, 

(3.45) Sd\t=o = Sd\ dn = 0. 
Let 

(3.46) 

G(t) = || I 



+ II S n u \\hp(a,t-M' 2 -i) + II d t S n u ||i,p(o,t;L9) 



+ II S n d || 2(1-1) + || S n d \\LP(0,t;W 2 >i) + II d t S n d \\LP(0,t:L") 

L~(0,t;B 5 , p p ) 

+ || V(5„P ||LP(0,t;i«) ■ 
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Repeating the argument from (3.22)-(3.37), we get 



(3.47) 



G(t)<-G{t) on [0,Ti]. 



Hence, G(t) = for t G [0,Ti], which implies the uniqueness on the interval 
[0,71]. 

To complete this section, we show the continuous dependence. Noticing the 
proof of Lemma 3.1 and Lemma 3.2, more precisely, from (3.19) and (3.36), we 
can deduce that for any given initial-boundary data, there exists a T\ > only 
depending on the initial-boundary data such that H(t) < 2C(1 + C 2 H )H n on 
[0,Ti]. Let 



(3.48) 
(3.49) 



(u ,d ) g d\J' p x (s 9 2 S Lp) nc 2 ' a (ffl)), 



2(1- 



(uJ,d5)->(uo,do) in D A , P ' x(B,; p p 'n^ a (9H)) 



Assume that (u n ,d n ) is the corresponding solution with the initial-boundary 
condition (itp , ) and («, d) is the solution with initial-boundary data (uo, do). 
Define 

(3.50) S n u = u-u n , S n d = d-d n , S n P = P-P n . 

Then, the triplet (S n u, S n d, S n P) satisfies the following system 



(3.51) 



dS n u 

at 



A5 n u + VS n P 



= -u ■ Vu + u n ■ Vu n - V(^^) 
+V(^f^) - Ad ■ Vd + Ad n ■ Ad n , 

= -u ■ Vd + u n Vd n 

+ I Vd I 2 d- I Vd n I 2 d n , 



V ■ S n u = 0, 

/ n S n Pdx = 

with initial-boundary conditions 



in (0,Ti) x 0, 

in (0,Ti) x 0, 
m (0,Ti) x O, 



(3.52) 
(3.53) 

Let 

(3.54) G n (t) 



5 n u\d P Q Tl =u - Uq, 
S n d\d P Q Tl = d - dg. 



II <HI 



L^(0,t-D A 



,p + II Su \\LP(0,t;W 2 -") + II d t 5u \\ L P(0,t;Li) 



+ || <W || 2(1 _i) + || Sd \\LP(0,t;W 2 -i) + II 9 t Sd \\ L p(0.t:Li) 

L~(0,t;S,, p p ) 

+ II V5P 



\LP(0,t;Li) ■ 

(3.55) G„(0) = \\uo-uS II i-i.p + ||do-d5 

^ AO 



2(1--) 



nc 2 <°(dQ) 
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Using Lemma 2.2 and Theorem 2.1, and repeating the proof from (3.29) to 
(3.36), we have that 

(3.56) G n (t)<CG n (0) + ^G n (t), on [O.Ti], 

where the constant C only depends on uq, do and T\. This implies G n (t) < 
CG n (0). Thus, we get the continuous dependence of the initial-boundary data. 
Now, we complete the proof of well-posedness. 



(4.1) 



4 Global existence 

In this section, our aim is to extend the local solution which is established in 
Section 3 to a global solution in case of small initial-boundary conditions. 

Given any unit vector e e S 2 such that do\d£l = e. Let (it, d) is a solution 
of the system (3.1)-(3.2). A directly computation shows that (u,d — e) satisfies 

^ + u ■ Vu - Ait + VP = -V • (Vd© Vd), in K+ x O, 

ff + it • Vd = (Ad+ | Vd | 2 d + e | Vd | 2 ), in R+ x fi, 

V • it = 0, in I + x!l, 
/ n Pda; = ) 

with initial-boundary conditions 

, . f (it(0,a;),d(0,x)) = (u (a:), d (a;) - e) x e 0, 

1 J \ (it(t, ar), x)) = (0, 0), (i, x) E M+ x 90. 

We first let T* < oo be the maximal time of existence for (u,d — e, P). Define 

(4.3) 

H{t) = || It || !_!,„ + || U \\ L v(0.t;W^.i) + || d t U \\ L P(0,t;Li) 

L°°(0.t:D Aq " ) 

+ l|d-e|| 2(1-1) + \\ 9 t (d- e) \\ L p/ ot . Lg \ + \\ P \\ LP / 0t . w i,g\ 

LP(0,t;B q , p p ) 
+ || d - e \\LP(0,t;W 2 -"), 

(4-4) 

Ha = II u o || i-i.p + || do - e || 2(1 _ i, . 

D v B v 

Aq a d,P 

Using the system (3.1), Theorem 2.1 and Lemma 2.2, we get 
(4.5) 

H(t) <C{H + || u • Vu ||iP ( o, ti L«) + II V • (V(d - e) V(d - e)) || 

+ || it • V(d-e) ||LP( ,t;L9) + HI V(d-e) | 2 (d- e) ||LP(o,t ; i«) + III V(d-e) | 2 e ||z, P (o,t ; z,«))- 
Noticing that | d | = | e |= 1 and (2.12)-(2.13), we have that 
(4.6) H(t) <C{H + || it || 2(1 _i) || u \\ L p(o,f,w^) 

+ II d - e || 2(i-i) | d - e \\Lp(o.t:W 2 .i) 

L°°(0,t;B q , p p ) 

+ II u II i-i, „ II d - e \\ L p(o.t.w 2 .i) 

+ II d - e || 2(i-i) | d - e \\Lp(o,t-,W 2 ."))- 
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From the definition of H{t) and (4.4), we have 

(4.7) H(t)<C{H +AH 2 (t)), on [0,T*). 

If we take Ho < 3 ^ ci , then the continuity and monotonicity of H (t) yields that 
there exists a positive T such that 

(4.8) H(t) < 2CH < j^, on [0,T]. 
On the other hand, (4.5) yields that 

»(.)< 1 -"-^' )i - «w> 1+ "- 8 f"- )i « io,n. 

Since i?(t) < j^yi on [0,T], the continuity of H(t) implies that 

( 4.io) *<,) < '-u-mc*.)* . 

This yields that 

(4-11) ||( u , d _ e ,P)|| B ,,.<^. 

It is impossible and the solution exists globally in time. This implies that 
(u, d, P) is a global solution of the system (3.1)-(3.2). 
The proof of Theorem 1.1. is complete. 
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